Dynamical properties of the synchronization transition 

Michel Droz 1 and Adam Lipowski 1 ' 2 

1 Department of Physics, University of Geneva, CH 1211 Geneva 4, Switzerland 
^Department of Physics, A. Mickiewicz University, 61-614 Poznan, Poland 

Abstract 

We study the dynamics of the synchronization transition (ST) of one-dimensional coupled map 
lattices. For the Bernoulli map it was recently found by Ahlers and Pikovsky (Phys. Rev. Lett. 88, 
254101 (2002)) that the ST belongs to the directed percolation (DP) universality class. Spreading 
dynamics confirms such an identification, only for a certain class of synchronized configurations. 
For homogeneous configurations spreading exponents w and 5 are different than DP exponents but 
their sum equals to the corresponding sum of DP exponents. Such a relation is typical to some 
models with infinitely many absorbing states. Moreover, we calculate the spreading exponents for 
the tent map for which the ST belongs to the bounded Kardar-Parisi-Zheng (BKPZ) universality 
class. Our estimation of spreading exponents are consistent with the hyperscaling relation. Finally, 
we examine the asymmetric tent map. For small asymmetry the ST remains of the BKPZ type. 
However, for large asymmetry a different critical behaviour appears with exponents being relatively 
close to the ones of DP. 
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I. INTRODUCTION 



Recently, synchronization of chaotic systems received considerable attention [|TJ. These 
studies are partially motivated by experimental realizations in lasers, electronic circuits, 
and chemical reactions 0. An interesting problem concerns synchronization in spatially 
extended systems ||. It turns out that in such systems synchronization can be regarded 
as a nonequilibrium phase transitions. The determination of the universality classes for 
nonequilibrium phase transitions is a problem much debated in the literature. Thus it is 
natural to ask whether the synchronization transition (ST) can be incorporated into an 
already known universality class. For certain cellular automata the synchronized state is 
actually an absorbing state of the dynamics, and ST in such a case was found to belong 
to the Directed Percolation (DP) universality class @j. This problem is more subtle for 
continuous chaotic systems as e.g., coupled map lattices (CML) ||. For CML a perfect 
synchronized state is never reached in finite time which weakens the analogy with DP. 
Recently Pikovsky and Kurths argued that synchronization of CML should generically be 
described by the so-called bounded Kardar-Parisi-Zheng (BKPZ) universality class || [7]]. 

However, it was observed that for noise-induced synchronization M as well as for a 
system of two interacting identical CML 0, the synchronization transition belongs to the 
DP universality class. It then was argued that DP critical behaviour might emerge for 
systems with discontinuous local maps, for which a linear approximation, that leads to the 
BKPZ universality class, breaks down ||. 

It is well known that models with a single absorbing state exhibit a DP criticality |IU |. At 
first sight one can think that the synchronized state is unique and that this is the reason of its 
relation with DP criticality in some cases. However, if one studies the dynamical properties 
of the ST by means of the so-called spreading approach [jllj , several synchronized states could 
be considered and some dynamical properties depend on the choice of a synchronized state. 
For example, the spreading exponents which we measure for homogeneous synchronized state 
take non-DP values. However, these exponents satisfy a more generalized scaling relation 
that also holds for some models with infinitely many absorbing states [|l!| . 

In addition to that, we measured spreading exponents for ST belonging to the BKPZ 
universality class. Our results, obtained for the symmetric tent map, together with previous 
estimations of other exponents in this universality class, show that the hyperscaling relation 
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is satisfied in this case. Finally, we examine the asymmetric tent map. For not too large 
asymmetricity the ST remains of BKPZ type. However, for large asymmetricity the nature 
of the ST changes and the critical exponents that we measured are relatively close to that 
of the DP. It shows that the non-BKPZ behaviour appears even for continuous maps. 

In section II we define the model and briefly describe the simulation method. The results 
are presented in section III and a final discussion is given in section IV. 

II. MODEL AND SIMULATION METHOD 

Our model is the same as the one examined recently by Ahlers and Pikovsky (AP) and 
consists of two coupled CML's |§, 

ui(x, t + 1)\ (l - 7 7 \ /(l + eA)f( Ul (x, t + 1)) 
u 2 (x, t + l)J ^7 1 - 7/ X \(1 + eA)/(w 2 (z, t + 1)) 
Av the discrete Laplacian Av(x) = v(x — 1) — 2v(x) + v(x + 1). Both space and time 

are discretized, x = 1, 2 ... ,L and t — 0, 1, Periodic boundary conditions are imposed 

ui,2(x + L,t) = ui t2 (x,t) and similarly to previous studies we set the intrachain coupling 
e = 1/3. Varying the interchain coupling 7 allows us to study the transition between 
synchronized and chaotic phases. Local dynamics is specified through a nonlinear function 
f(u) and several cases will be discussed below. 

We introduce a synchronization error w(x, t) = \ui(x, t) — ii2{x, t)\ and its spatial average 
w(t) = r J2t=i w ( x , t). The time average of w(t) in the steady state will be simply denoted 
as w. In the chaotic phase, realized for sufficiently small 7 < 7c , one has w > 0, while in 
the synchronized phase (7 > 7c ), w = 0. Moreover, at criticality, i.e. for 7 = j c , w(t) is 
expected to have a power-law decay to zero w(t) ~ t~ e . 

It is well known that the spreading dynamics is a very effective method to study phase 
transitions in models with absorbing states fll], pL3| . In this method one prepares the model 
in an absorbing state and then one locally sets the activity and monitor its subsequent 
evolution. Typical observables in this method are the number of active sites N(t), the 
probability P(t) that activity survives at least up to time t and the average square spread 
of active sites R 2 (t). However, as we already mentioned, our model never reach a perfectly 
synchronized state in a finite time. To define P(t) we have to introduce a small threshold p 
and consider a system as "perfectly" synchronized when w(t) < p. As we will show below 

3 




estimation of the critical value of the interchain coupling 7 C as well as some exponents seems 
to be independent on the precise value of p, as long as it remains small. 

As for the other observables, we define them in the following, ^-independent way 



where xq denotes the site where the activity was set initially. 

According to general scaling arguments |nj, we expect that at criticality (i.e., for 7 = 7 C ): 
w(t) ~ t v , P(t) ~ t~ s , and R 2 (t) ~ t z . These relations define the critical exponents r], 5, 
and z. 

Note that in general the critical behavior of w(t) ~ t v in the spreading case might be 
different from the non-spreading one w(t) ~ t~ e discussed above. 

III. RESULTS 

A. Bernoulli map 

First, let us consider the case when f(u) is a Bernoulli map, namely f(u) = 2w(mod 1) 
and < u < 1. Measuring w in the steady state and its time dependence w(t), AP concluded 
that the synchronization transition in this case takes place at 7 = 7 C = 0.2875(1) and belongs 
to the DP universality class. 

Our purpose here is to apply the spreading dynamics to this map. First, we have to set 
the model in a synchronized state i.e., the variables of both chains of CML's must take the 
same value. But despite that constraint there is a considerable freedom in doing that. Below 
we present results for two choices: 

(i) Random synchronized state 

In this case each pair of local variables takes a different random value, i.e., Ui(x, 0) = 
U2(x, 0) = r(x), where < r(x) < 1 and r(x) varies from site to site. 

(ii) Homogeneous synchronized state 

In this case r{x) = r and is constant for each x. 

Having set the model in the synchronized state we initiate an activity assigning at a randomly 
chosen site Xq a new random number either for Ui(xq) or for ^2(^0)- Then we monitor the 
time evolution, measuring w(t). We always used the system size L = 2t m + 1 where t m is 





4 




FIG. 1: The survival probability P(t) for the Bernoulli map at 7 = 7 C = 0.2875 and for random 
and homogeneous absorbing states as a function of time t. Calculations were done for the threshold 
p = 10 -10 (dashed line) and 10~ 16 (solid line). The upper straight dotted line have a slope which 
corresponds to the DP value 5 = 0.1595 and for the lower line 5 = 0.43. 

the maximum simulation time. In such a way the spreading is never affected by finite size 
effects. Moreover, the data were averaged over 10 3 — 10 4 independent runs. Final results 
shown in Fig. |] - Fig. || are obtained for 7 = 7 C = 0.2875. Off-criticality the data deviate 
from the power law behaviour in a typical way. Estimating the asymptotic slope of our data 
we conclude that for the random synchronized state the exponents are in a good agreement 
with very accurately known DP values: 5dp = 0.1595, rjop = 0.3137, and z = 1.2652 fll4"| . 



However, in the case of homogeneous synchronized states, the exponents rj and S are clearly 
different and we estimate 5 = 0.43(1) and 77 = 0.05(1). 

Let us notice that a very similar situation occurs in certain models with infinitely many 
absorbing states, where these critical exponents are also found to depend on the choice of 
an absorbing state []P2"| . Although the spreading exponents are nonuniversal in this class of 
models, they obey the following relation 



77 + 5 = r] DP + S DP . 



(3) 



It turns out that the estimated exponents for homogeneous synchronized states also obey 
this relation within an estimated error. From Eq. (§) and the hyperscaling relation [13| 



v + s + e = ^ (d = i) 



(4) 



it follows that the exponent z must be constant and thus independent on the choice of the 




FIG. 2: The averaged difference w(t) as a function of time t for the Bernoulli map at 7 = j c = 
0.2875 and for random and homogeneous absorbing states. Calculations were done for the threshold 
p = 10 -10 (dashed line) and 10 -16 (solid line). The upper straight dotted line have a slope which 
corresponds to the DP value 77 = 0.3134 and for the lower straight line r/ = 0.05. 



absorbing state. Indeed, in Fig. ^ one can see that both for the random and homogeneous 
absorbing states the asymptotic increase of R 2 (t) is described by the same exponent. 

The main conclusion which follows from the above calculations is that synchronization 
should be considered as a transition in a model with multiple absorbing states rather than 
a single one. Except for non-DP values of rj and 5 it has essentially no other consequences 
in the one-dimensional case. For example all steady-state exponents keep the DP values. 
However, the situation might be different in higher-dimensional models with infinitely many 
absorbing states. In particular, there are some analytical and numerical indications that in 
such a case a non-DP criticality might appear (TH . 



B. Symmetric Tent map 

The second map examined by AP is the symmetric tent map denned as f(u) = 1 — 2\u — 
1/2 1, (0 < u < 1). In this case they found that the synchronization transition belongs to 
the BKPZ universality class and occurs at 7 = 7 C = 0.17614(1). A qualitative difference 
between the critical behaviour in the case of Bernoulli and tent maps is attributed to the 
strong nonlinearity (i.e., discontinuity) of the first map. Another interesting feature reported 
by AP is a different behaviour of the transverse Lyapunov exponent Aj_. For the tent map 
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FIG. 3: The average spread square R 2 (t) as a function of time t for the Bernoulli map at 
7 = 7 C = 0.2875 with random (solid line) and homogeneous (dashed line) absorbing states. The 
straight dotted line have a slope which corresponds to the DP value z = 1.2652. 



Aj_ vanishes exactly at j c , while for the Bernoulli map it vanishes inside a chaotic phase. 

For the tent map we also performed the spreading-dynamics calculations but only for 
random synchronized states. Our results for the time dependence of w(t) are shown in 
Fig. |. The asymptotic decay of w(t) is described by the exponent r\ = —0.53(5) (note a 
minus sign). 

As for the estimation of 5 our results strongly suggest that in this case 6 = 0. Indeed, in 
the examined range of the threshold p we observed that all runs survived until the maximum 
measured time t m = 10 4 . From the behaviour of R 2 {t) (Fig. [5]) we estimate z = 1.30(5), 



which is in a good agreement with the expected KPZ value z 



EM 



The hyperscaling relation provides a constraint on the exponents rj, 6, O, and z. Since for 
the BKPZ universality class z is a simple number (4/3), it is tempting to assume that other 
exponents are also simple numbers. Thus we might speculate that in this model the exact 
value of r\ is -0.5. Then, from the hyperscaling relation (^) together with 6 = and z = | 



we obtain = ^ = 1.166 . . .. Such a value might be compared with numerical estimations 
of this exponent which range from 1.1(1) [0, [HJ to 1.26(3) 
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FIG. 4: The averaged difference w(t) as a function of time t for the symmetric tent map at 
7 = 7 C = 0.17614 and for random absorbing states. Calculations were done for the threshold 
1CT 10 . The least-square fit to the last decade date give 77 = —0.53(5). The dotted straight line has 
a slope which corresponds to 77 = —0.5. 




log 10 (t) 

FIG. 5: The average spread square R 2 (t) as a function of time t for the symmetric tent map 
at 7 = 7 C = 0.17614 and random absorbing states. The straight dotted line have a slope which 
corresponds to the KPZ value z = 1.3333. 

C. Asymmetric Tent map 

In this subsection we examine the asymmetric tent map denned as 

{au for < u < 1/a 

(5) 
a(l - u)/(a - 1) for 1/a < u < 1, 

and 1 < a < 2. For a = 2 this is the symmetric tent map. Let us notice that in the limit 
a — > 1, the slope of the second part of this map diverges. Below we present the results of 
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FIG. 6: The steady-state synchronization error w as a function of 7 for the asymmetric tent map 
(a = 1.1). The measurement was made during t = 10 5 steps after 5 • 10 4 steps of relaxation. In 
the range (0.22 < 7 < 7 C = 0.2288) w is well fitted with a power-law function a{^ c — 7) , where 
= 1.5. 

our calculations for a = 1.1, 1.02 and 1.01. 

1. a=l.l 

When a is not too close to 1, model (P remains in the BKPZ universality class. In 
particular, for a = 1.1 the estimated exponents (3 = 1.5(1) (Fig. |6]) and 9 = 1.16(10) 
(Fig. [7]) are in good agreement with other estimations for this universality class. 

2. a=1.02 



However, for a close to 1 the nature of the transition changes. From the measurements of 
w (Fig. H) and its time dependence w(t) (Fig. |^) we estimate that in this case 7 C = 0.17095 
and O = 0.25(5). Using the spreading dynamics calculations for the random synchronized 
states at the critical point we estimate (Fig. [10] — Fig. [12|): 77 = 0.15(3), 5 = 0.27(2), and 
z = 1.32(3). All of the exponents, except z are far from the BKPZ values. 
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FIG. 7: The time dependence of w for the asymmetric tent map (a = 1.1) and 7 equal to (from 
top) 0.228, 0.2285, 0.2288 (critical), and 0.229. Calculations were done for L = 3 • 10 5 . The dotted 
straight line has a slope corresponding to = 1.16. 
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FIG. 8: The steady-state synchronization error w as a function of 7 for the asymmetric tent map 
(a = 1.02). The measurement was made during t = 10 5 steps after 5 • 10 4 steps of relaxation. 

3. a=1.01 



From the measurements of w (Fig. O) and its time dependence w(t) (Fig. Ill) we estimate 
that in this case 7 C = 0.13695 and G = 0.15(2). Similarly to the a = 1.02 our results for w 
close to the critical point are not sufficiently accurate to determine (3. Using the spreading 
dynamics calculations for the random synchronized states at the critical point we estimate 
(Fig. 0-Fig. [T|): 77 = 0.22(2), 5 = 0.21(2), and z = 1.25(5). Again, the exponents, except 
z are far from the BKPZ values. However, G and z are consistent with the DP value. 
Moreover, the sum 77 + 5 = 0.43 that is also quite close to the DP value (0.4732). Thus, it is 
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FIG. 9: The time dependence of w for the asymmetric tent map (a = 1.02) and 7 equal to (from 
top) 0.1707, 0.1709, 0.17095, 0.1710, and 0.1712. Calculations were done for L = 5 • 10 4 . We 
identify the central curve as critical. 
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FIG. 10: The averaged difference w(t) as a function of time t for the asymmetric tent map and 
for random absorbing states. Calculations were done for a = 1.02 with 7 = 7 C = 0.17095 and for 
a = 1.01 with 7 = 7 C = 0.13695. The dotted straight lines have a slope which corresponds to 
T) = 0.15 (upper) and r] = 0.22 (lower). 

in our opinion likely that in this case the model belongs to the DP universality class but the 
random synchronized states are not "natural" absorbing states |12[] and that is why r] and S 
separately do not take their DP values. Let us notice that for the Bernoulli map, examined 
in section III-A, the random synchronized states are probably a good approximation of the 
"natural" absorbing states, since the obtained spreading exponents coincide the DP values. 
Our numerical results are summarized in Table |. 
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FIG. 11: The survival probability P(t) for the asymmetric tent map with random absorbing states. 
Calculations were done for a = 1.02 with 7 = 7 C = 0.17095 and for a = 1.01 with 7 = j c = 0.13695. 
The dotted straight lines have a slope which corresponds to 6 = 0.27 (upper) and 6 = 0.21 (lower). 
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FIG. 12: The average spread square R 2 (t) as a function of time t for the asymmetric tent map and 
random absorbing states. Calculations were done for a = 1.02, 7 = 7 C = 0.17095 (solid line) and 
a = 1.01, 7 = 7c = 0.13695 (dashed line). The straight dotted lines have slopes that correspond 
to DP and KPZ value of z. 

IV. DISCUSSION 

In this paper we studied the dynamical properties of CML models undergoing a synchro- 
nization transition. Our main results show that:(i) the synchronized state is not unique 
and (ii) non-BKPZ critical behaviour might appear for continuous maps. One of the open 
questions that we leave for the future is why for sufficiently asymmetric but continuous 
maps the linear analysis, that leads to the relation with BKPZ model, breaks down. For 
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FIG. 13: The steady-state synchronization error w as a function of 7 for the asymmetric tent map 
(a = 1.01). The measurement was made during t = 10 5 steps after 5 • 10 4 steps of relaxation. 
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FIG. 14: The time dependence of w(t) for the asymmetric tent map (a = 1.01) and 7 equal to 
(from top) 0.1368, 0.13685, 0.1369, 0.13695, and 0.137. Calculations were done for L = 5 • 10 4 . 



TABLE I: Critical parameters for the symmetric and asymmetric tent map. The last column shows 
how much the critcal exponents deviate from critical hyperscaling relation (|3|). 
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7] + 5 


§-e-77-<5 


2 


0.17614 


1.26(3) 


-0.53(5) 


0.0(1) 


1.50(5) 


1.30(5) 




0.07 (0.11) 


1.1 


0.2288 


1.16(10) 






1.5(1) 


1.30(5) 






1.02 


0.17095 


0.25(5) 


0.15(3) 


0.27(2) 
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1.32(3) 


0.42 


0.01 (0.11) 


1.01 


0.13695 


0.15(2) 


0.22(2) 


0.21(2) 


? 


1.25(5) 


0.43 


0.04 (0.08) 


DP 




0.1595 


0.3137 


0.1595 


0.2765 


1.2652 


0.4732 
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strongly asymmetric maps the (deterministic) noise is very large and this seems to be the 
only factor that might change the universality class. It is known that depending on how 
the noise scales with the order parameter, Langevin-type models, that presumably describe 
our model at criticality and at a coarsed-grained level might exhibit either DP or BKPZ 
critical behaviour |T7[. In particular when the amplitude of noise scales linearly with the 



order parameter the model exhibits the BKPZ critical behaviour while the DP one for the 
square root scaling. Thus, the problem is to explain how a qualitatively different scaling 
of the noise in the Langevin-type model could emerge due to a more quantitative changes 
(in asymmetricity) in a coupled CML system. It is also tempting to expect that BKPZ and 
DP critical behaviours meet for a certain asymmetricity value a c (1.1 < a c < 1.01) at a new 
multi-critical point. Our estimation of O for a = 1.02 significantly deviates both from DP 
and BKPZ values which might be an indication of such a new critical point. 
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